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Abstract 

The tau neutrino with a mass of about 10 MeV can be the "late decaying particle" 
in the cold dark matter scenario for the formation of structure in the Universe. We 
show how this may be realized specifically in the recently proposed doublet Majoron 
model. 



Study of the mass and interaction of neutrinos is a long-standing subject in particle 
physics and offers one of the key clues to possible new phenomena beyond the standard 
model. Majoron models have been attracting a lot of interest in this respect. They provide 
the neutrinos with Majorana masses and a new type of interaction not present in the standard 
model. The interaction is due to the coupling of the neutrinos and other matter to the 
Majoron, a Nambu-Goldstone boson associated with the violation of the lepton number 
symmetry. This extra interaction was originally used to facilitate the decay of a massive 
neutrino which would otherwise be ruled out due to the cosmological constraint on the 
neutrino mass [0. In a different application, it allows a stable massive neutrino to be 
suitable as a dark matter candidate 0. In this paper, we examine another possibility that 
this interaction may bring about; the r neutrino (z/t-) as a candidate for the late decaying 
particle in the cold dark matter (CDM) scenario for the formation of structure in the Universe 




The idea of the late decaying particle was proposed to reconcile a setback of the CDM 
model in explaining the formation of large-scale structure in the Universe |^ ^. This setback 
became more evident by the recent COBE detection of anisotropy in the temperature of the 
cosmic background radiation: The theoretical prediction for an f2 = 1 inflationary Universe 
on the power spectrum of the density fluctuation gives a larger power than the observation 
at small scales A < lOh^^ Mpc (the Hubble constant Hq = lOOh km/s/Mpc ) once it is 
normalized at large scales A ~ 10^h~^ Mpc using the COBE detection P|. A remedy is to 
delay the time of matter-radiation equality, at which sub-horizon-sized fluctuations begin to 
grow. The delay reduces the power on the small scales, while the resulting larger horizon 
size means relatively more large-scale structure than the simple CDM scenario. 

A massive particle species that decays into relativistic particles can do this trick. Since 
the decreasing rate of the relic energy density of massive matter is smaller than that of 
radiation, its energy necessarily dominates the Universe if it is sufficiently long-lived. Its 
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subsequent decay into relativistic particles gets the Universe into the radiation-dominant era 
again with more energy, and delays the time of forthcoming matter-radiation equality. 

In the late decaying particle scenario, we have two matter-radiation equality epochs: first 
the relic dominates the Universe and then the cold dark matter does. They are separated 
by a radiation-dominated era after the z/T--decay. We distinguish values of cosmological 
variables at these epochs with subscripts 'eqi' for the former and 'eq2' for the latter, e.g., 
the age of the Universe tEQi, the temperature Teqi, the horizon size Aeqi, etc. 

Recently, Dodelson et al. examined a scenario in which a massive r neutrino, with its 
mass in a range m^^ ~ 1 - 10 MeV, may be the late decaying particle [Q. (See also Ref. 
for other candidates in particle physics models.) An important constraint for this scenario is 
the one from primordial nucleosynthesis: the equivalent number of massless neutrino species 

should be less than 3.3 or 3.04 but a r neutrino of this mass range may possibly 



contribute more than this bound 10 



This difficulty is rather easily evaded in a Majoron model, thanks to the new interaction 
between v-r and the Majoron ipi^. In the standard model, the relic abundance Yoo of J^r, the 



ratio of its number density to the entropy density (see |n| for the definition), is determined 
by the speed of the annihilation process via Z° exchange compared to the Universe expansion 
||Tl| , [12| . In Majoron models the process 

v^v^ ^ (pi^ipi^, (1) 

also works to decrease the relic abundance. The coupling is proportional to m^^/vi,, 

where is a scale for the lepton number violation. Thus the process (|l|) can be still active 
even after the exchange process shuts off if this m,^^ to ratio is sufficiently large. Then 
the relic density can be much smaller than the one in the standard model. This can make 
invisible with respect to the dynamical evolution of the Universe at the time of primordial 
nucleosynthesis and avoid the resulting constraint. In other words, the lifetime upper bound 



of 100 seconds estimated for a heavy to be the late decaying particle in Ref. is no 
longer a constraint in this case. 

The doublet Majoron (DM) model which we have proposed recently is very suitable 
for the late decaying particle scenario. An advantage of this model, compared with the 
singlet Majoron model, is that it allows us to use a smaller scale fL for the lepton number 
violation and thus a larger z/^ - coupling without conflicting with the basic concept of 
the seesaw mechanism We estimate Teqi for various values of vi^ and my^, and show 
that it indeed provides a possibility that z/^ can be the late decaying particle. 

Let us first estimate how low the temperature Teqi should be in order to satisfy the 
constraint from primordial nucleosynthesis. The precise definition for Teqi is the tempera- 
ture at which the r neutrino energy p^,^ = (27r^/45)mj,^Foo fl'*s T^ becomes the same as the 
radiation energy = {tt^ / 30) g^T^, 

Teqi = — m^^Yoo, 2 

45 5f* 

where g^ and g^s are the statistical weights of the light degrees of freedom for the energy 
and entropy density, respectively. They are functions of temperature; specific values of g^ 
and g^:S at a given temperature depend on details of the thermal history of the Universe. 
The values we specifically use correspond to the following situation: Ur was so heavy, m,^^ 
> a few MeV, that it was non-relativistic when the Z° exchange process for its annihilation 
shuts off; its abundance is frozen at the temperature Tf, which is assumed to take place 
before electrons (e) and positrons (e"*") annihilate in pairs, i.e., Tf > a few tenths of MeV. 
Thus we use g^ = g^s = 10 (a sum of contributions from photon (7), e, e"*", z/e, u^, and ipi^) 
for evaluating Tf and Foo- After the pair annihilation of e and e"*", they reduce to (7* = 3.17 
and g^^s = 3.64. If z/^ is relativistic at the termination of the Z° process, these values change. 
Also there may occur a deviation in the temperatures of 7 and ipi^. We, however, neglect 
these subtleties in this paper. For lighter neutrino, m^^ ~ 1 MeV, this may cause an error; 
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but the numerical error in the final results of Teqi expected by this neglect is at most 30 - 
40 % and does not affect the discussions we make in this paper. 

The temperature at which primordial nucleosynthesis commences is about 1 MeV [ |11| |. 
At this temperature, we need to satisfy the condition (7^ < 11.3 or 10.82, which correspond 
to the bounds < 3.3 |^ or 3.04 respectively. Since = 10 at T ~ 1 MeV, the Vr 
fraction in the energy density must satisfy 

^ < 0.08. (3) 

PR 

On the other hand, the fraction in terms of Teqi is given by 

Pv^ g*si^ MeV) fi'*(rEQi) Teqi Teqi 



(4) 



PR g,{l MeV) ^.s(Teqi) 1 MeV 1 MeV 

as long as Tf is higher than 1 MeV. Thus a r neutrino with Teqi less than 10~^ MeV is safe 
in this respect. 

To evaluate Fqo, we use formulas given by Kolb and Turner 



^ 3.79(n + l)xr^ 

and 

Xi = ln[0.038(?2 + l)(5f/5f^,^/^)mpim,,^ao] 

-{n + 1/2) ln{ln[0.038(r2 + l){g/g.y^)mpim,^ao]}, (6) 

where mpi = 1.2 x 10^^ GeV is the Planck mass, g = 2 for z/^, n and cxo are read off from the 
form of the thermally averaged cross section for the process 

{cTV,,,)=aoiT/m,X- (7) 
The freezing temperature is given by Tf ~ m^^xY^- 
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The evaluation of the cross section is straightforward. The appropriate interaction terms 
are given by 



L = — (z/^V^z/^) H — (cos a 0_ + sin a I 



COS a 0_ + sin a { ( z/t- ia Vr) — [ J^r icr v^* ) } 



2t;L 

where are neutral scalar fields, a and (3 are mixing angles, which are the same as those 
defined in Ref. The corresponding cross section in the singlet Majoron model has been 



evaluated in Refs. [|T^. We calculate the same Feynman diagrams as those in [|1^: two 



diagrams with the Vt - (fi^ vertices and two scalar-exchange diagrams. The freezing of the 
abundance mostly takes place when the initial two z/^'s are non-relativistic and their energy 
is much smaller than the mass of (j)±. In this energy region, the 0±-exchange diagrams are 
negligible and we obtain the same result as the one in Ref. 



a=— (9) 
967r E ^ ^ 

in terms of the energy E and the momentum p of one of the initial neutrinos in the center 
of mass frame. A characteristic behavior of a, i.e. it vanishes at vanishing \p\, comes 
from a combined effect of the statistics of identical particles and the conservation of both 
angular momentum and "CP" ; the latter can be defined as a discrete symmetry of the Ur - 
ipi^ interaction and forbids the s-wave contribution. The thermally averaged cross section is 
given by the integration over the distribution function for the relative velocity frei = 2|p 



This result gives ctq = (l/327r)(my^^/t>L'^) and n = 1 in (0). 

^ Note that the initial neutrinos in the process (|l|) are identical Majorana particles. Thus the event rate 
of the process per unit comoving volume is (1/2) of (ctWioi) defined this way. In the Boltzman equation, 
which is the basis to derive the formulas (||)-(||) this factor 1/2 is cancelled by another factor 2 that 
represents the fact that two neutrinos annihilate in the process. 
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We plot the values of Teqi obtained by Eqs. (Q), (|), and (||) for various values of 
and m^^ as contours in Fig. 1. The thermal history of the Universe we have assumed based 
on the DM model is correct for most of the range of values of rriu^ and fL shown in Fig. 1. 
Obviously Teqi needs to be higher than Teq2, which is about 1 eV, in order that u^- plays the 
role of the late decaying particle. Thus the parameter region in Fig. 1 for accommodating 
the late decaying particle is 

(Note that z/^ with Teqi less than 1 eV can be a candidate of dark matter [Q].) 

In the DM model, smaller than 10 GeV predicts a light scalar boson 0_, which can be 
a rare decay product in Z*^ — 0_ + (a fermion pair), and contradicts the known lower mass 
bound, about 60 GeV, for the standard Higgs boson |[T3|lp|- Thus a region fL > 10 GeV is 
left for the r neutrino to be the late decaying particle. The temperature Teqi is, then, 1 - 
10 keV. 

The lifetime of the r neutrino, t^^, should be adjusted so that it generates an appropriate 
amount of radiation energy in its decay. Let us estimate the required lifetime. We use 
the sudden-decay approximation and assume Uj- decays all at once at the age — 
and temperate Td. The total radiation energy density after the decay, including the decay 
product R', is 

To fit the predicted power spectrum of the density fiuctuation to the observation, this needs 
to be about 3 times bigger than the radiation energy in the standard prediction |p. Thus 

Teqi/Td — 2. Taking into account the relation T oc t"^/^ in the i/T--dominated era, we get 

^Since ^_-Z'^-Z°coupling is proportional to (— cos /3 sin a + sin /3 cos a), we can fix this problem by tuning 
the two mixing parameters in the model, i.e., a ~ but we will not pursue it in this paper. 



7 



Tj,^ ~ StEQi- Since the age at Teqi is 



tEQi ^ 2.4 X 10^9 ( ) sec, (13) 
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^EQl 



where Tq = 2.735i^' is the present (photon) temperature of the Universe, r^^is 10'' - 10^ 
seconds. The hfetime in the DM model is given by ||l3| 



647r 

where Ry^u^ is the flavor changing matrix element between Vj. and lighter neutrinos (z/^ = 
fe, J^/i)- If we take — 20 GeV and rriu^ ~ 10 MeV to get an idea of the magnitude of \R\, 
it resides in a range \R\ ~ 10^^ - 10^^'^. The smallness of these values can be regarded as 
a result of the seesaw mechanism: \R\ can be parametrized as [{my^/M) sin^] with M the 
mass scale of the gauge singlet neutrino and 9 a mixing angle; the values we used for rriy^ 
and imply M ~ 10 TeV and \R\ ~ 10~*^ 9. 

In a late-decaying-particle scenario for structure formation, we necessarily have an extra 
small scale corresponding to the horizon at tEQi- Its size Aeqi at Teqi, after being scaled up 
to the present taking into account the expansion of the Universe, is 

Aeqi ^ 4.8 x 10^ (^] Mpc. (15) 

V^EQl/ 

Thus Aeqi — 10 - 100 kpc for the allowed parameters in the DM model. The consequence 
of the existence of this scale for structures in the Universe needs to be clarifed by further 
investigations. It may be related to the dwarf galaxies P]. 

In summary, we have shown that the DM model represents an appropriate particle physics 
model for realizing a possibility that the r neutrino, with a mass of about 10 MeV and a 
lifetime in the range 10^ to 10^ seconds, is the late decaying particle in the CDM scenario 
for the formation of structure in the Universe. 
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Figure caption 

Fig. 1 Contour plot of log[TEQi/lGeV]. The contours correspond to log[rEQi/lGeV] = 
—5, —6, —7, —8, —9, —10 from top to bottom. 



10 
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